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ABSTRACT 

Dust in the host galaxies of gamma-ray bursts (GRBs) dims and reddens their 
afterglow spectra. Knowledge of the nature of this dust is crucial for correcting 
for extinction, providing clues to the nature of GRB progenitors, and probing 
the interstellar medium of high-redshift galaxies as well as the nature of cosmic 
dust when the universe was much younger and galaxies were much less evolved. 
The dust and extinction properties of GRB host galaxies are still poorly known. 
Unlike previous work, we derive in this Letter the extinction curves for 10 GRB 
host galaxies without a priori assumption of any specific extinction types (such as 
that of the Milky Way, or the Small/Large Magellanic Clouds). It is found that 
there appears to exist two different types of extinction curves: one is relatively 
flat and gray, the other displays a steeper dependence on inverse wavelength, 
closely resembling that of the Milky Way but with the 2175 A feature removed. 

Subject headings: dust, extinction — gamma rays: bursts 

1. Introduction 

The existence of dust extinction towards gamma-ray bursts (GRBs) and their afterglows 
has been well established observationally through (1) "dark bursts" - nearly ~60% of the X- 
ray afterglows have no optical counterparts - the undetected optical afterglows may have been 
extinguished by dust in the host galaxy (see Lazzati et al. 2002 and references therein); (2) 
reddening - the optical/near-infrared (IR) spectral energy distributions (SEDs) of afterglows 
deviate from that expected from standard models, 1 indicative of dust reddening (see Kann 
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-'As any host-galaxy extinction would result in a steepening of the intrinsic spectral slope, a steeper 
(observed) slope (than all intrinsic spectral slopes) betrays the presence of host galaxy dust. 
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et al. 2006 and references therein); (3) the reduced Balmer line ratios (as a consequence of 
dust extinction) in the spectra of some GRB host galaxies (e.g. see Djorgovski et al. 1998) 
compared with the expected ratios for the standard Case B recombination (Osterbrock & 
Ferland 2006); (4) the depletion of dust-forming heavy elements such as Si and Fe in some 
host galaxies (e.g. see Savaglio et al. 2003); and theoretically (5) through the association of 
(at least long-duration [>2s]) GRBs with massive stars and star-forming regions embedded 
in dense clouds of dust and gas (Paczyhski 1998). 

A precise knowledge of the nature (e.g. size, composition) of the dust in GRB host 
galaxies is very useful for (1) correcting for the extinction of afterglows from X-ray to near- 
IR wavelengths; (2) constraining the nature of the GRB progenitors (i.e. collapsing massive 
stars or merging neutron binaries); and (3) understanding the interstellar medium (ISM) of 
high-redshift galaxies and the cosmic star formation history (e.g. see Ramirez-Ruiz et al. 
2002). 

However, our knowledge of the dust in GRB host galaxies is very limited. Previous 
studies all are limited to a qualitative and rather speculative analysis of the extinction curves 
derived by fitting the observed optical/near-IR SEDs of GRB afterglows with a power- law 
(approximating their intrinsic spectra) reddened by an extinction law of known types (mostly 
either that of the Milky Way [MW], the Small Magellanic Cloud [SMC], or that of the Large 
Magellanic Cloud [LMC]) adopted as a priori (e.g. see Stratta et al. 2004, Kann et al. 2006). 

In this Letter, by assuming that the GRB afterglows have a power-law intrinsic spectrum 
as expected from the standard fireball model, we obtain the extinction curves for 10 GRB 
host galaxies without a priori assumption of the extinction law. The size and composition of 
the dust are modeled quantitatively in terms of the silicate-graphite interstellar grain model. 

2. Method 

In the standard fireball model, the temporal (t) and frequency (v) dependence of GRB 
afterglows can be well described by F„(t) oc t~ c V~ /3 , where both a and f3 are related with 
the electron energy distribution index p (dn oc E~ p dE for the energy distribution of the 
shock-heated electrons; see Sari et al. 1998). Therefore, with the decay index a determined, 
p as well as the intrinsic spectral index (5 can be well constrained. 

In addition to the extinction of the host galaxy, GRB afterglows are also subject to the 
extinction of our own MW galaxy. We correct for the latter using the reddening maps of 
Schlegel et al. (1998). The host's spectrum is redshifted to longer wavelengths where any 
errors in the Schlegel et al. map are not important. 
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Let F° be the intrinsic spectrum. With the Galactic reddening corrected, the ob- 
served wavelength (A)-dependent spectrum is F\ = F° exp (— A\), where A\ is the ex- 
tinction from the host galaxy. Setting V band as the zero point, we obtain Fx/ Fy = 
(Ay/A) 2 "* exp (A v - A x ), where A and Ay = 5500 A are both measured in the rest frame. 
With (3q estimated from the decay index a, we can obtain A\ — Ay from Fx/Fy. To de- 
termine Ay, we fit a simple function A x — Ay — a/ [l + (A/A^) 2 ] + b for A > Ay and then 
extrapolate this function to A — > oo to get Ay (——b) since A x — > as A — >■ oo. 

The standard silicate-graphite interstellar dust model, consisting of a mixture of spher- 
ical silicate and graphite grains, is shown successful in reproducing the extinction and IR 
emission of the MW galaxy, SMC and LMC (Weingartner & Draine 2001; Li & Draine 2001, 
2002). We will apply this grain-mixture to model the extinction curves Ax/Ay determined 
for the GRB host galaxies, but with a simpler functional formula for the dust size distri- 
bution: dn ~ a~ v exp(— a/a c ) da for both grain types, where a is the grain radius, ranging 
from 0.05 to 2.5 /xm, n is the power-law index, a c is the cut-off size, and dn is the number of 
grains in the size interval [a, a + da}. The mass fraction of graphite dust is / gra [for silicates 
it is (1 — /gra)]- m Figure 1 we show that the MW, SMC and LMC extinction curves are 
well fit by the silicate-graphite model with this simple dust size distribution function. 

3. Data 

To construct an extinction curve, one requires simultaneous multi-band photometry. 
We carefully collected such data from literature for ten bursts taken when their decays 
were in a steady power-law state(e.g. see Panaitescu & Kumar 2001, Fan & Piran 2006 for 
very detailed analysis). The optical and near-IR UBVRIJHK fluxes F v of these bursts are 
tabulated in Table 1. Also tabulated is flo, the intrinsic spectrum index calculated from 
fitting the decay index based on the standard model of afterglow (see §2). 

To compare the extinction properties in the rest frame of each burst, in the observer 
frame we need to calculate Ax — A\yt 1+Z \]. The results are shown in Table 2. We correct the 
wavelength with a factor (1 + z) to get Ax — Ay in the rest frame. Also shown in Table 2 
are Ay, obtained by extrapolating Ax — Ay to A — > 00 (see §2). 

4. Extinction Curves: Two Different Types? 

The extinction curves (normalized to the V band) derived for these 10 host galaxies are 
shown in Figures 2 and 3. A common feature for these curves is the absence of the 2175 A 
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bump, which is the strongest absorption band in the MW extinction (Li 2005). They appear 
to fall into two categories: one is relatively gray, showing a much flatter increase with inverse 
wavelength A -1 (we call it "Type-I"; see Fig. 2), the other displays a much steeper, almost 
linear increase with A^ 1 (we call it "Type-II"; see Fig. 3). 

In view of its possible practical use (e.g. in correcting for extinction of afterglows), 
following Fitzpatrick & Massa (1990), we take an analytical fit to these 2 extinction types: 

/ ci + c 2 x + c 3j D(x;7,x ) + c 4 F(x), £>l/xm _1 , 
Ax,Av -\kx^\ x<l^-\ (1) 



where x = \ l , D(x; 7, xq) = x 2 / (x 2 — x^) 2 + x 2/ y 



, and 



0.5392 (x - 5.9) 2 + 0.05644 (x - 5.9) 3 , x > 5.9 /^m" 1 , 
[X) ~\0, x<5.9fim-\ {2} 

In this fit, Ci and Ci determine a linear "background" term; C3 determines the strength of the 
2175 A feature which is represented by the "Drude profile" term D(x; 7, x ) - the theoretical 
profile for a classic damped harmonic oscillator (Bohren & Huffman 1983); 7 and Xq are 
respectively the FWHM and peak position of the Drude profile; 2 C4 determines the far-UV 
curvature term represented by F(x). The resulting analytical fits are plotted in Figures 2 
and 3; and the fitted parameters c±, C2, C3, C4, and 7 are tabulated in Table 3. 

To quantify the dust properties, we fit both extinction types in terms of the silicate- 
graphite dust model. The model fits and parameters are respectively shown in Figures 2,3 
and in Table 4. The model for the flat "Type-I" extinction curve is dominated by large 
silicate grains (with / gra « and a flat size distribution). In contrast, the model for the 
steep "Type-II" extinction has a much steeper size distribution, indicating the richness of 
small dust in the host galaxies with a "Type-II" extinction curve. 



5. Discussion 

Previous efforts in deriving the dust extinction of GRB host galaxies all assume an initial 
power-law for the GRB intrinsic spectrum (as expected from the standard fireball model), 
and an extinction curve of known types such as that of MW, SMC, LMC, the Calzetti et al. 
(1994) law suitable for local starburst galaxies, and the Maiolino et al. (2001) law suggested 



2 In the Galactic ISM, the strength and width of the 2175 A extinction bump vary with environment while 
its peak position is quite invariant (see Li 2005). We therefore fix x at 4.6 ^imT 1 . 
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for the dust in the circumnuclear region of AGNs (e.g. see Stratta et al. 2004 and Kann et 
al. 2006 for recent examples). Our approach, without the need for an initial assumption of 
specified extinction types, is more favourable because of the lack of a priori knowledge of the 
extinction in the GRB hosts. The fact that neither of the two extinction types derived in this 
work resembles that of MW, SMC and LMC (see Figs. 2,3) challenges the initial assumption 
of specific extinction types. 

In consistent with previous studies, this work also finds no evidence for the 2175 A 
extinction bump. 3 The carrier of this bump remains unidentified over 40 years after its 
first detection. It is generally believed to be caused by aromatic carbonaceous (graphitic) 
materials, likely a cosmic mixture of polycyclic aromatic hydrocarbon (PAH) molecules (see 
Li 2005 and references therein). The nondetection of the 2175 A extinction bump in the 
GRB hosts could be explained in terms of (1) the depletion of its carrier (e.g. PAHs) by 
condensation onto the ice mantles coated on the dust in the dense star-forming regions in 
which GRBs are embedded, and (2) the destruction of its carrier by UV/X-ray radiation in 
the immediate vicinities (up to ~20pc) of bursts. 

Our Type-I extinction curve (see Fig. 2) is more gray (i.e. with a much weaker depen- 
dence on wavelength) than that of MW, SMC, and LMC. Gray dust has been invoked by 
a number of authors (e.g. Savaglio et al. 2003, Savaglio & Fall 2004; Stratta et al. 2004, 
2005) to account for the small ratios of extinction and/or reddening (derived from fitting 
the afterglow SEDs) to H column densities (determined from X-ray or Lya absorption) of 
GRB hosts, which are usually smaller than that of MW by a factor of ~ 10-100 (Galama & 
Wijers 2001, Simon et al. 2001, Hjorth et al. 2003, Vreeswijk et al. 2004). 4 Gray dust could 
be created by (1) the preferential destruction of small grains by the intense UV and X-ray 
radiation from the GRB (Waxman & Draine 2000, Fruchter et al. 2001, Perna et al. 2003) in 
the immediate GRB environment, say, within ~ 10-20 pc around the burster (Hjorth et al. 
2003, Savaglio et al. 2003, Stratta et al. 2005), and (2) the growth of dust through coagula- 
tion (which also leads to the depletion of small grains) in the high-density environments such 



3 The only exception is GRB 991216. Its afterglow spectrum shows a depression in flux in between 4000 A 
and 5500 A, suggesting the presence of a red 2175 A- type extinction bump (at ~ 2360 A) in the host of GRB 
991216 (Vreeswijk et al. 2006). Similar red 2175 A bumps have been previously seen in a few UV-strong, 
hydrogen-poor stars in the MW galaxy (see Li 2005 and references therein). 

4 Alternatively, these results can be explained by an intrinsicaly low dust content - a low metallicity and/or 
a low level of depletion of the dust- forming heavy elements in the burst environment (Fynbo et al. 2003, 
Hjorth et al. 2003, Vreeswijk et al. 2004, Fruchter et al. 2006). However, Savaglio et al. (2003) found that in 
the host galaxies of 3 GRBs both the column densities of metals (as indicated by Zn) and the depiction of 
heavy elements (such as Fe, Si and Cr) are large, indicating a large dust content. 
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as the cores of star- forming regions harboring GRB events (Maiolino et al. 2001, Stratta et 
al. 2004). Both mechanisms naturally lead to a dust size distribution skewed toward large 
grains, producing an extinction curve weakly dependent on wavelength. 

Nongray extinction has also been reported (e.g. for the host of GRB 010222 by Galama 
et al. 2003). Our nongray, Type-II extinction curve (see Fig. 3) shows a steep increase with 
A -1 , indicating that small grains are more abundant in the GRB host galaxies with a Type- 
II extinction than in those with a Type-I extinction. This is probably caused by a less 
complete destruction of small grains (e.g. associated with less energetic bursts of less intense 
UV/X-ray radiation), and/or a slower coagulational growth of dust (e.g. associated with less 
dense environments or metal-poorer galaxies). Since our sample is small, we do not want to 
overinterpret its implications. A more systematic investigation of the extinction and dust 
properties of a large sample of GRB host galaxies is in progress. 

Finally, if one really wants to make a priori assumption of the extinction spectral shape 
when deriving the GRB host galaxy dust extinction, we suggest the adoption of the Fitz- 
patrick k Massa (1990) formulae instead of that of MW, SMC or LMC. By first visually 
inspecting the GRB spectrum, one can set c 3 = (in eq.[l]) if there is no sign of the 2175 A 
extinction bump or curvature on large wavelength intervals. 
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comments. S.L.C. and D.M.W. are supported by the NSFC grants 10225314 and 10233010, 
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Fig. 1. — Fitting the extinction curves of the Milky Way galaxy, the SMC and LMC by the 
silicate-graphite model with a simple dust size distribution dn/da oc a~ v exp(— a/a c ). Points 
- observational data; solid lines - model fits. 
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Fig. 2. — Flat, "Type-I" extinction curves (normalized to the V band) for 4 GRB host 
galaxies. Both axes are in the rest-frame of the GRB host. Also plotted are the model 
(black) and analytical (magenta) fits to the observed extinction curves. The MW, SMC and 
LMC extinction curves are also shown for comparison. 
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Fig. 3. — Same as Figure 2 but for the relatively steep, "Type-IP extinction curves derived 
for 6 GRB host galaxies. 
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Note. — Data are all taken from the cited references. The Galactic visual extinction w is derived from 
the reddening maps of Schlegel et al. (1998). The UBVRIJHK fluxes are measured in the observer frame; 
z is the redshift of the burst; t is the time when the data was taken since the burst; and fio is the intrinsic 
spectral index derived from the standard afterglow model (see §2). 

References. — (1) Frail et al. 2003; (2) Jensen et al. 2001; (3) Jakobsson et al. 2003; (4) Stratta et al. 
2005; (5) Covino et al. 2003; (6) Holland et al. 2003; (7) Klose et al. 2004; (8) Vreeswijk et al. 2004; (9) 
Bloom et al. 2004; (10) Jakobsson et al. 2004. 
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Table 2: Extinction in the rest frame of each burst calculated from Table 1. Ay is obtained 

by extrapolating an analytical fit of A\ — Aiv(i+ Z )] to A — > oo (see §2). 

GRB z J A v A x - A mi+Z)] 

U B V R I J H K 

980703 0.966 L07 (148 ~ 029 029 CL24 4L41 -0.26 

000301C 2.0404 0.56 0.35 0.51 0.43 0.34 0.26 0.21 0.1 0.01 -0.13 

011211 2.14 0.63 0.07 0.18 0.12 0.13 0.07 0.09 -0.002 ••• 0.002 

020405 0.691 0.5 0.57 0.66 0.47 0.30 0.28 0.04 -0.39 -0.45 -0.51 

020813 1.255 0.52 0.43 0.65 0.52 0.48 0.39 0.31 -0.01 -0.26 -0.09 

021004 2.33 0.45 0.51 1.42 0.99 0.65 0.55 0.43 ••• 0.06 

030226 1.986 0.33 0.40 ••• 0.45 0.36 0.28 0.22 0.09 -0.002 -0.20 

030323 3.3718 0.28 0.92 ••• 2.64 1.71 1.29 0.99 0.84 0.64 0.15 

030329 0.1685 0.5 0.36 0.14 0.08 -0.05 0.03 -0.16 -0.26 

030429 2.658 0.63 0.41 0.74 0.52 0.42 0.14 ••• -0.03 
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Table 3: Parameters for fitting the extinction curves of GRB host galaxies with the Fitz- 
patrick & Massa (1990) formulae. Note that x is fixed at 4.6 pm" 1 and is not a free 
parameter (see Eq.[l]). 
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Table 4: Parameters for modeling the 2-type extinction curves of GRB host galaxies (as well 
as that of MW, SMC and LMC) with the silicate-graphite dust model. 



Galaxy r] a c (//m) / gra X 2 /dof 
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a X 2 = EaiiAEaiiGRBs K A \ / A v) mod - { A x/ A v ) oh f I <? 2 is obtained by summing up all wavebands and all 
GRBs, where a is the uncertainty for a given GRB at a given band. 



